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Abstract 

This paper is devoted to the prediction problem in extreme value theory. Our 
main result is an explicit expression of the regular conditional distribution of a max- 
stable (or max-infinitely divisible) process {r]{t)}t£T given observations = 
Hi, 1 < i < k}. Our starting point is the point process representation of max- 
infinitely divisible processes by Gine, Hahn and Vatan (1990). We carefully analyze 
the structure of the underlying point process, introduce the notions of extremal 
function, sub-extremal function and hitting scenario associated to the constraints 
and derive the associated distributions. This allows us to explicit the conditional 
distribution as a mixture over all hitting scenarios compatible with the conditioning 
constraints. This formula extends a recent result by Wang and Stoev (2011) dealing 
with the case of spectrally discrete max-stable random fields. This paper offers new 
tools and perspective for prediction in extreme value theory together with numerous 
potential applications. 
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1 Introduction 
1.1 Motivations 

Since the pioneer works by Fisher and Typett [11] and Gnedenko piS], the univariate 
theory of extremes is now well established with extensive studies on models, domains 
of attraction, parameter estimations, etc. (see e.g. de Haan and Fereira [8] and the 
references therein). The last decades have seen the quick development of multivariate and 
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spatial extreme value theory: the emphasis is put on the characterization, modeling and 
estimation of the dependence structure of multivariate extremes. Among many others, 
the reader should refer to the excellent monographs [21 [HI [ini [E] and the reference therein. 

Max-stable random fields turn out to be fundamental models for spatial extremes since 
they arise as the the limit of rescaled maxima. More precisely, consider the component- 
wise maxima 

rjnit) = max t eT, 

l<i<n 

of independent and identically distributed (i.i.d.) random fields {Xi(t)}t^T, i > I. If the 
random field rjn = {f]n(t)}teT converges in distribution, as n — )■ oo, under suitable affine 
normalization, then its limit rj = {r]{t)}t^T is necessarily max-stable (see e.g. [HI HI])- 
Therefore, max-stable random fields play a central role in extreme value theory, just like 
Gaussian random fields do in the classical statistical theory based on the Central Limit 
Theorem. 

In this framework, the prediction problem arises as an important and long-standing 
challenge in extreme value theory. Suppose that we already have a suitable max-stable 
model for the dependence structure of a random field rj = {?7(t)}teT and that the field is 
observed at some locations ti, . . . ,tk G T. How can we take benefit from these observa- 
tions and predict the random field t] at other locations ? We are naturally lead to consider 
the conditional distribution of {ri{t)}t^T given the observations {f]{ti) = Vi, 1 < < k}. 
A formal definition of the notion of regular conditional distribution is deferred to the 
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In the classical Gaussian framework, i.e., if 77 is a Gaussian random field, it is well 
known that the corresponding conditional distribution remains Gaussian and simple for- 
mulas give the conditional mean and covariance structure. This theory is strongly linked 
with the theory of Hilbert spaces: the conditional expectation, for example, can be ob- 
tained as the L^-projection of the random field rj onto a suitable Gaussian subspace. In 
extreme value theory, the prediction problem turns out to be difficult. A first approach 
by Davis and Resnick [3 IB] is based on a L^-metric between max-stable variables and 
on a kind of projection onto max-stable spaces. To some extent, this work mimics the 
corresponding L^-theory for Gaussian spaces. However, unlike the Gaussian case, there 
is no clear relationship between the predictor obtained by projection onto the max-stable 
space generated by the variables {rjiti), 1 < i < k} and the conditional distributions of 
r] with respect to these variables. A first major contribution to the conditional distri- 
bution problem is the work by Wang and Stoev p5j. The authors consider max-linear 
random fields, a special class of max-stable random fields with discrete spectral measure, 
and give an exact expression of the conditional distributions as well as efficient algo- 
rithms. The max-linear structure plays an essential role in their work and provides major 
simplifications since in this case r] admits the simple representation 

r]{t) = \/ Z,f,{t), teT, 
i=i 

where the symbol \/ denotes the maximum, fi,---,fq are deterministic functions and 
Zi, . . . , Zq are i.i.d. random variables with unit Frechet distribution. The authors deter- 
mine the conditional distributions of (-^j)i<j<g given observations {rjiti) = yi, I <i < k}. 
Their result relies on the important notion of hitting scenario defined as the subset of 
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indices j G ll,q} such that rjiti) = Zjf(ti) for some i G |1, k], where, for n > 1, we note 
[Ij-n,] = {1, . . . ,n}. The conditional distribution of {Zj)i<j<q is expressed as a mixture 
over all admissible hitting scenarios with minimal rank. 

The purpose of the present paper is to propose a general theoretical framework for 
conditional distributions in extreme value theory, covering not only the whole class of 
sample continuous max-stable random fields but also the class of sample continuous max- 
infinitely divisible (max-i.d.) random fields (see Balkema and Resnick [Ij). Our starting 
point is the general representation by Gine, Hahn and Vatan [12] of max-i.d. sample 
continuous random fields (see also de Haan [7J for the max-stable case). It is possible to 
construct a Poisson random measure $ = Xlili ^<t>i the space of continuous functions 
on T such that 

N 

vit) = \JUt), teT. 

Here the random variable N is equal to the total mass of $ that may be finite or infinite 
and = stands for equality of probability laws (see Theorem [l] below for a precise state- 
ment). We denote by [$] = {(^j, 1 < i < N} the set of atoms of $. Clearly, 0(t) < rjit) 
for all t G T and G [$]. The observations {ri{ti) = i/i, I < i < k} naturally lead to 
consider extremal points: a function G [$] is called extremal if (j){ti) = ri{ti) for some 
i G otherwise it is called sub-extremal. We show that under some mild condition, 

one can define a random partition B = {6i, . . . ,6e) of {ti, . . . ,tk} and extremal func- 
tions ipf, . . . ,(p^ G [$] such that the point ti belongs to the component 9j if and only if 
^'jiti) = Tjiti). Using the terminology of Wang and Stoev [16], we call hitting scenario a 
partition of {ti, ■ ■ ■ , t^} that reflects the way how the extremal functions yj^, . . . , hit 
the constraints v^^(tj) < ''^(tj), 1 <i <k. The main results of this paper are Theorems 
andjsj where the conditional distribution of r] given {//(t?) = Ui, 1 < i < fc} is expressed 
as a mixture over all possible hitting scenarios. 

The paper is structured as follows. In Section 2, the distribution of extremal and sub- 
extremal functions is analyzed and a characterization of the hitting scenario distribution 
is given. In Section 3, we focus on conditional distributions: we compute the conditional 
distribution of the hitting scenario and extremal functions and then derive the conditional 
distribution of rj. Section 4 is devoted to examples: we specify our results in the simple 
case of a single conditioning point and consider max-stable models. The proofs are 
collected in Section 5 and some technical details are postponed to an appendix. 

1.2 Preliminary on max-i.d. processes 

Let T be a compact metric space and C = C(T, M) be the space of continuous functions 
on T endowed with the sup norm 

11/11 = sup |/(t)|, /gC. 

Let (fi, J-", P) be a probability space. A random process rj = {ri(t)}t^T is said to be max- 
i.d. on C if r/ has a version with continuous sample path and if, for each n > 1, there 
exist {rjni, I < i < n} i.i.d. sample continuous random fields on T such that 

n 

^ = V 

1=1 
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where Y denotes pointwise maximum. 

Gine, Hahn and Vatan (see [I2j Theorem 2.4) give a representation of such processes 
in terms of Poisson random measure. For any function f on T and set A C T, we note 
/(A) = sup,g^/(t). 

Theorem 1. (Gine, Hahn and Vatan JT3]) 

Let h be the vertex function of a sample continuous max-i.d. process rj defined by 

hit) = snp{x e R; P(r/(t) > x) = 1} E [-oo, oo), t E T, 

and define Ch = {f E C; f ^ h, f > h} . Under the condition that the vertex function 
h is continuous, there exists a locally-finite Borel measure fi on Ch, such that if ^ is a 
Poisson random measure $ on Ch with intensity measure fi, then 

{^WW = {sup{/^(t),0(t);0G (1) 

where [$] denotes the set of atoms of ^. 
Furthermore, the following relations hold: 

h{K) = sup{x E R; F{ri{K) > x) = 1}, K cT closed, (2) 

and 

P [r]iK,) <x„ l<i<n]= exp[-/i {U'l^.if E /(iT,) > x,})], (3) 
where n E N, Ki G T closed and Xi > h{Ki), 1 < i < n. 

Theorem [T] provides an almost complete description of max-i.d. continuous random 
processes, the only restriction being the continuity of the vertex function. Clearly, the 
distribution of r] is completely characterized by the vertex function h and the so called 
exponent measure fi. The random process e"^ — is continuous and max-i.d. and its 
vertex function is identically equal to 0. Since the conditional distribution of rj is easily 
deduced from that of — e^, we can assume without loss of generality that /i = 0; the 
corresponding set Co is the space of non negative and non null continuous functions on 
T. 

We need some more notations from point process theory (see Daley and Vere- Jones 
[21 Hj). It will be convenient to introduce a measurable enumeration of the atoms of $ 
(see [1] Lemma 9.1.XIII). The total mass of $ is noted = <l>(Co). If /x(Co) < oo, has 
a Poisson distribution with mean //(Co), otherwise = +oo almost surely (a.s.). One 
can construct Co- valued random variables (0i)i>i such that $ = Yld=i 

Let Mp(Co) be the space of point measures M = Ylii^^i^fi such that 

{fi E Co: II /ill > e} is finite for all e > 0. 

We endow Mp(Co) with the cr-algebra M.p generated by the applications 

M t-^ M(A), A C Co Borel set . 

For M E Mp(Co), let [M] = {fi, i E 1} he the countable set of atoms of M. If M is 
non null, then for all t E T, the set {f{t);f E [M]} is non empty and has finitely many 
points in (e, -|-oo) for all e > so that the maximum max{/(t);/ E [M]} is reached. 
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Figure 1: A representation of the point process $ (left) and of the associated maximum 
process 1] = max($) (right) in the moving maximum max-stable model based on the 
gaussian density function. Here T = [0, 5]. 

Furthermore by considering restrictions of the measure M to sets {/ G Cq; ||/|| > e} and 
using uniform convergence, it is easy to show that the mapping 



is continuous with the convention that max(M) = if M = 0. 

In Theorem [l] (with h = 0), Equation ([s]) implies that the exponent measure fi satisfies, 
for all e > 0, 



point measure $ and of the corresponding maximum process 1] = max($) in the moving 
maximum max-stable model based on the Gaussian density function. 

2 Extremal points and related distributions 

In the sequel, t] denotes a sample continuous max-i.d. random process with vertex func- 
tion h = and exponent measure fi on Cq. On the same probability space, we suppose 
that a Mp(Co)-valued Poisson random measure $ = YliLi ^4>i with intensity measure /i is 
given and such that rj = max($). 

2.1 Extremal and sub-extremal point measures 

Let C T be a closed subset of T. We introduce here the notion of i^'-extremal points 
that will play a key role in this work. We use the following notations: if /i, /2 are two 
functions defined (at least) on K, we write 



Let M e Mp(Co). An atom / G [M] is called i^'-sub-extremal if and only if / <x 
max(M) and i^-extremal otherwise. In words, a sub-extremal atom has no contribution 
to the maximum max(M) on K. 





fi =K /2 if and only if 
fi <K f2 if and only if 
/i /2 if and only if 



Vt G K, /i(t) = Ut), 
Vt G K, h{t) < Ut), 
3t G K, f\{t) > Mt). 
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Figure 2: Decomposition of the Poisson point measure $ into the i^'-extremal point 
measure $^ (black) and the i^'-sub-extremal point measure (grey). Left: K = [0,5]. 
Right: K = {3} represented by a black square. 

Definition 1. Define the K-extremal random point measure $^ and the K- sub -extremal 
random point measure $^ by 

N N 
i=l i=l 

Figure provides an illustration of the definition. It should be noted that and 



$^ are well defined measurable random point measures (see Lemma |3j in Appendix A.3). 
Furthermore, it is straightforward from the definition that 

$ = <I>J + max($;^) =K V and max($;^) <k V- 



Define the following measurable subsets of Mp(Co) (see Lemma |4] in Appendix A.3): 

= {m G Mp(Co); V/ G [M], / max(M)}, (5) 

CM = {MeMp{Coy,Wfe[M],f<Kg], (6) 

where g is any continuous function defined (at least) on K. Clearly, it always holds 

$J G and <l>^ G Cj^{r]). 

The following theorem characterizes the joint distribution of ($X' ^k) given that $^ 
is finite. We note So the Dirac mass at 0. 

Theorem 2. For all measurable A,Bg Mp(Co), 

eAxB, <l>+(Co) = 0] = exp[-/x(Co)]5o(^)5o(5), 

and, for k > 1, 

^[{<l>i,<^K)eAxB, $+(Co) = A;] 

Theorem [2] fully characterizes the joint distribution of ($X' ^k) provided that $^-(Co) 
is almost surely finite. We now focus on this last condition. 
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Proposition 1. 

The K -extremal point measure $J is a.s. finite if and only if one of the following condi- 
tion holds: 

(i) /i(Co) < +00; 

(ii) /u(Co) = +00 and inf ri(t) > almost surely. 

It should be noted that any simple max-stable random field (with unit Frechet mar- 
gins) satisfies condition (ii) above. See for example Corollary 3.4 in |12j . 

Remark 1. Using Theorem^ it is easy to show that the distribution of ($^, $^) has 
the following structure. Define the tail functional fix by 

Uk{9) = ^({/ e Co; / 9}) 

for any continuous function g defined (at least) on K. Suppose that $J is finite almost 
surely. Its distribution is then given by the so-called Janossy measures (see e.g. Daley 
and Vere- Jones section 5.3). The Janossy measure of order k of the K- extremal point 
measure $^ is given by 

Jk{dfu ...,dfk)= exp [-ftK (Vti/i)] gc+1/^®''(^/i' • • • '^/fc)- 

Furthermore, given that $^ = X]^=i'^/,:j conditional distribution of is equal to 
the distribution of a Poisson random measure with measure intensity lf^j^^yk_^j.yfi{df) . 
These results are not used in the sequel and we omit their proof for the sake of brevity. 

2.2 Extremal functions 

Let t E T. We denote by /it the measure on (0, +00) defined by 

fXt{A) = /i({/ e Co; fit) e A}), A C (0, +00) Borel set, 
and by fit the associated tail function defined by 

fit{x) = fit{[x,+ 00)), x>0. 

Note that 

P(77(t) < x) = exp(-/i({/ G Co; /(t) > x})) = exp(-/^t(x)), x > 0. (7) 

The following proposition states that, under a natural condition, there is almost surely a 
unique {t}-extremal point in $. This extremal point will be referred to as the t-extremal 
function and noted 0^. 

Proposition 2. For t E T, the following statements are equivalent: 

(i) $|j|(Co) = 1 almost surely; 

(ii) /it(0+) = +00 and fit is continuous on (0, +00); 
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(iii) the distribution of r]{t) has no atom. 



If these conditions are met, we define the t- extremal function (pf by the relation $^ 



a.s.. For all measurable A d we have 



P(0+gA) = / exp[-/i,(/(t))]/i(rf/). 

J A 



An important class of processes satisfying the conditions of Proposition [2] is the class 
of max-stable processes (see section 4.2 below). 



2.3 Hitting scenarios 

Proposition [2] gives the distribution of $J when K = {t} is reduced to a single point. 
Going a step further, we consider the case when K is finite. In the sequel, we suppose 
that the following assumption is satisfied: 

(A) K = {ti, . . . , tk} is finite and, 
for all t E K, fit is continuous and /if(0"'") = +00. 

Roughly speaking, this ensures that the maximum r]{t) = max($)(t) is uniquely reached 
for all t & K. This will provide combinatorial simplifications. More precisely, under 
Assumption (A), the event 

t&K 

is of probability 1 and the extremal functions 0^ , . . . , J are well defined. In the next 
definition, we introduce the notion of hitting scenario that reflects the way how these 
extremal functions hit the maximum rj on K. 

Let Vk be the set of partitions of K. It is convenient to think about K as an ordered 
set, say ti < ■ ■ ■ < t^. Then each partition r can be written uniquely in the standardized 
form r = {ti,...,Ti) where i = £(r) is the length of the partition, Ti G K is the 
component of ti, r2 C -ft' is the component containing mm{K \ Ti) and so on. With this 
convention, the components ri, . . . , of the partition are labeled so that 

min Ti < ■ ■ ■ < min T£. 

Definition 2. Suppose that Assumption (A) is met. Define ~ the (random) equiva- 
lence relation on K = {ti, . . . ,tk} by t t' if and only if (f)f = (p'f,. The partition 
= (6'i, 6*^(0)) of K into equivalence classes is called the hitting scenario. For 
j G |1,£(B)], let (fj be the extremal function associated to the component 6j, i.e., such 
that (f^ = (f)f for all t G 6j. 

We illustrate the definition with two examples in Figure [s] Clearly a point G [$] is 
i^-extremal if and only if it is t-extremal for some t G -ft', so that = {0^, t G K}. 
Furthermore, the random measure is almost surely simple, i.e. any atoms have a 
simple multiplicity, otherwise the condition $|"^^(Co) = 1 a.s. would not be satisfied for 
some t G K. These considerations entail that 

£(6) 

i=i 
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Figure 3: Two realisations of the Poisson point measure $ and of the corresponding 
hitting scenario and extremal functions yj^, . . . , ^9^^^ with K = {ti, ^2, ^3, ^4} repre- 
sented by the black squares. Left: the hitting scenario is 9 = ({^i}, {^2}, {^3, ^4}), the 
extremal functions are ipf = (f)f_^, = 4>ti and '^t — 4>ti = (f^ti- Right: the hitting 
scenario is 9 = ({ti,t2}, {^35^4}), the extremal functions (black) are '^t = (pfi = (Pt2 
¥^2 = <Pt 



In particular, the length i{Q) of the hitting scenario is equal to $^(Co). Furthermore 
the extremal functions satisfy 

Vj G Vt e 9„ ^^(t) > V v^+(t). (10) 

The distribution of the hitting scenario and extremal functions is given by the following 
proposition. The proof relies on Theorem [2] 

Proposition 3. Suppose Assumption (A) is met. 

Then, for any partition r = (ri, . . . , r^) G Vk, and any Borel sets A C Cg, B C Mp(Co); 
we have 



= r, ((^+,...,(^+) G A ^-^eB] (11) 

= / h^mm, v,v, e 5 n Q (vj^ j,)] ^^''\dh^ • • • , ^^)- 

J A 

3 Regular conditional distribution of max-id processes 

We now focus on conditional distributions. We will need some notations. 

If s = (si, ...,si) G and / G Cq, we note /(s) = (/(si), . . . , f{si)). Let be the 
exponent measure of the max-i.d. random vector ?7(s), i.e. the measure on [0, +00)' \ {0} 
defined by 

/is(A) = ^({/ G Co; /(s) G A]), A C [0, +CX))' \ {0} Borel set. 
Define the corresponding tail function 

/is(x) = /i({/ G Co; /(s) ^ x}), x G [0, +00)' X {0}. 
Let {-Ps(x, df); x G [0, +00)' \ {0}} be a regular version of the conditional measure /i((i/) 



given /(s) = x (see Lemma [2] in Appendix A.2). Then for any measurable function 

F : [0,+oo)' X Co ^ [0,+cx)) 
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vanishing on {0} x Cq, we have 

/ F(/(s),/)/i(rf/) = / / F(x,/)P,(x,t;/K(rfx). (12) 

^Co i[0,+oo)'\{0} JCo 

Let t = {ti, . . . and y = {yi, . . . ,yk) G [0, +00)^^. Before considering the condi- 
tional distribution of rj with respect to //(t) = y, we give in the next theorem an exphcit 
expression of the distribution of ?7(t). We note K = {ti, . . . ,1^}. For any non empty 
L C -ft', we define L = {i e |1, A;] : U e L} and set ti = {U)i(zL, Yl = {yi)ieL and 
L^ = K\L. 

Theorem 3. Suppose assumption (A) is satisfied. For t G Vk, define the measure 
on [0, +00)^^ by 

iyl{C) = P(r/(t) G C; e = r), C C [0, +00)'= Borel set. 

Then, 

I 

ulidy) = exp[-/it(y)] (g) {Pt.^. (y.,, {/(t..) < y.c}) /i^.^ (rfy.j} (13) 

i=i 

and the distribution of ri{t) is equal to Vt = TIit^^Vk 

Under some extra regularity assumptions, one can even get an explicit density function 



for Vt (see the section 4.3 on regular models below). 

We are now ready to state our main result. In Theorem |4] below, we consider the 
regular conditional distribution of the point process $ with respect to ?7(t) = y. Then, 
thanks to the relation rj = max($), we deduce easily in Corollary [s] below the regular 
conditional distribution of 77 with respect to ?7(t) = y. 

Recall that the point process has been decomposed into two parts: a hitting scenario 
together with extremal functions {ipf, . . . , V^^g^) and a i^'-sub-extremal point process 
Taking this decomposition into account, we introduce the following regular conditional 
distributions: 

7rt(y, ■) = P[0G-h(t)=y] 
gt(y,r, •) = P[(^+)G-h(t)=y,e = r] 
^t(y,r, (/,),-) = P[<|.^G-h(t)=y,e = r,(y,;) = (/,)]. 

We use here the short notations i = ii^r), (y?^) = {{pf , . . . , ip'^) and similarly (/j) = 
{fi, . . . , fe). The following theorem provides explicit expressions for these regular condi- 
tional distributions. 

Theorem 4. Suppose assumption (A) is satisfied. 
1. For any t G Vk, it holds Vt{dy)-a.e. 

vrt(y,r) = ^(y) (14) 

where ut and vl are defined in Theorem^ and dul/dut denotes the Radon-Nykodym 
derivative of w.r.t. Ut. 
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2. It holds h't{dy)nt{y,dT)-a.e. 



(15) 



/n words, conditionally on 'q{t) = y anc? 6 = r, t/ie extremal functions {ipf, . . . , y^^) 
are independent and ipj' follows the distribution P^^ (y^-^,, df) conditioned to the con- 
straint fi^r'r) < Yt't- 

3. Let C-(y) = {M G Mp(Co); V/ G [M], /(t) < y}. It holds a.e. 

«.(y...(/,).fl)-i^.(y.s)^^l|J|^ (16) 



for any measurable B G Mp(Co). /n words, conditionally on rj(i) = y, $^ zs 
independent of G anc? (v^i",- • • ,V^^e)) <J?^<i /^os the same distribution as a Poisson 
point measure with intensity l{j(t)<y}/i((i/). 

As a consequence, we deduce the regular conditional distribution of rj with respect to 

^(t) = y- 

Theorem 5. It holds vt{dy)-a.e. 

P[r/(s) < z I r^it) = y] = exp[-/i({/(s) ^ z, /(t) < y})] 

^Pt./y.,,{/(t.;)<y.;, /(s)<z}) 

^ .S- n ^.^(y,,U(t,)<y,}) 

/or any / > 1, s G T' and z G [0, +oo)'. 

Remark 2. Let us mention that Theorem |4] suggests a three-step procedure for sampling 
from the conditional distribution of rj given ri{t) = y: 

1. Draw a random partition r with distribution '7rt(y, ■). 

2. Given r = {ri, . . . , r^}, draw i independent functions ipi, . . . ,'ipi, with ipj following 
the distribution PtrXyTi,df) conditioned on /(t^c) < y^c. 

J •' 3 3 

3. Independently of the above two steps, draw Ylii^i ^ Poisson point measure on Cq 
with intensity l{/(t)<y}/i((i/). It can be obtained from a Poisson point measure with 
intensity fi{df) by removing those points not satisfying the constraint /(t) < y. 

Then, the random field 

fi(t) = max{^pi{t), . . .,i^i(t)} V max{0j(t), i e I}, t eT, 

has the required conditional distribution. 

The issues and computational aspects of conditional sampling are addressed in the 
paper [9]. The special case of Brown- Resnick max-stable processes is considered and 
tractable expressions are derived and the above three-step procedure is implemented 
effectively. 
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4 Examples 



As an illustration, we apply in this section our general results to specific cases. 

4.1 The case of a single conditioning point 

It is worth noting that the case of a single conditioning point, i.e. k = 1, gives rise to 
major simplifications. There exists indeed a unique partition of the set K = {t} so that 
the notion of hitting scenario is irrelevant. Furthermore, there is a.s. a single i^-extremal 
function (pf which is equal to the t-extremal function 0^. In this case. Theorems [4] and 
[5] simplify into the following proposition. 

Proposition 4. Lett E T and suppose that conditions (i)-(iii) in Proposition^ are met. 
Then, conditionally on ri{t) = y, (pf and are independent; the conditional distribution 
of ipf is equal to Pt{y,-); the conditional distribution of ^^^y is equal to the distribution 
of a Poisson point measure with intensity Furthermore, for / > 1, s G 

and z G [0, +oo)', 

P[r;(s) < z I r^{t) = y] 
= P,{y, {/(s) < z}) exp[-/i({/(s) ^ z, f{t) < y})]. (17) 

4.2 Max-stable models 

We put the emphasis here on max-stable random fields. For convenience and without 
loss of generality, we focus on simple max-stable random fields rj, i.e., with standard unit 
Frechet margins 

P(r/(t) <x)= exp[-a;-^]l{,.>o}, x G M, t G T. 
A random field t] is said to be simple max-stable if for any n > 1, 

n 
i=l 

where {r^j, i > 1} are i.i.d. copies of r]. Any general max-stable random field can 
be related to such a simple max-stable random field r] by simple transformation of the 
margins, see e.g. Corollary 3.6 in [12]. Furthermore, Corollary 4.5.6 in [8j states that rj 
can be represented as 

m}t^T = {\jTm)] ^ (18) 

i>i *^ 

where (rj)j>i is the nonincreasing enumeration of the points of a Poisson point process on 
(0, oo) with intensity x~'^dx, (Fj)j>i is an i.i.d. sequence of continuous random processes 
on T, independent of (Ti)i>i and such that 

E[Yi{t)] = l, teT, and E[||Fi||] < oo. 



12 



Since a continuous simple max-stable random field is max-i.d., it has a Poisson point 
measure representation ([T]). The normalization to unit Frechet margins entails that the 
vertex function h is equal to and that the exponent measure /i satifies, for all t eT, 

Mdy) = y''^^{y>o}dy and fitiv) = y'\ y >o. 

The correspondence between the two representations ([T]) and (18) is the following: the 
point measure $ = ^j>i Sr^Yi is a Poisson point measure on Co with intensity 

POO 

/i(yl) = / E[l{ry^g^}] r~'^dr, A C Co Borel set, 
Jo 

The distribution of the Fj's, denoted by a, is called the spectral measure and is related 
to the exponent measure /i by the relation 

poo r 

= / / l^rfeA} (j{df)r^'^dr, A C Co Borel set. 
Jo Jco 

Taking into account this particular form of the exponent measure, we can relate the 
kernel Pt{y, ■) to the spectral measure a. For x G M, we note (x)"*" = max(x,0). 

Proposition 5. Let t] be a continuous simple max-stable random field with spectral mea- 
sure a and t & The {t}- extremal function (pf has conditional distribution 

P[0+ G ■ I ii{t) = y] = Ptiy, •) = ^ l{^^f^.yf{t)a{df). 

Furthermore, for / > 1, s G and z G [0, +ooy , 

F[r]{s) < z I r]{t) = y] (19) 

expT f /(^^) /Wa^ 



^ ^(rf/) / 1 n^^mf{t)a{df). 



Equation (19) extends Lemma 3.4 in Weintraub [T7] where only the bivariate case 
Z = 1 is considered. Note the author considers min-stability rather than max-stability; 
the correspondence is straightforward since, if rj is simple max-stable, then ri~^ is min- 
stable with exponential margins. 



4.3 Regular models 

We have considered so far the case of a single conditioning point which allows for major 
simplifications. In the general case, there are several conditioning points and the hitting 
scenario is non trivial. This introduces more complexity since the conditional distribution 
is expressed as a mixture over any possible hitting scenarios and involves an abstract 
Radon-Nykodym derivative. The framework of regular models can be helpful to get more 
tractable formulas. 

The exponent measure is said to be regular (with respect to the Lebesgue measure) 
if for any / > 1 and s G with pairwise distinct components, the measure fj^s{dz) is 
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absolutely continuous with respect to the Lebesgue measure c/z on [0, +oo)^ We denote 
by hs the corresponding Radon-Nykodym derivative, i.e., yUs('^z) = hs{z)dz. 

Under this assumption, we can reformulate Theorems [3] and |4| For example. Equation 
(13) implies that the distribution z/t of //(t) is absolutely continuous with respect to the 
Lebesgue measure with density 



dy 



(y) = exp[-/it(y)] 5^ n / hitr^,t,c){yr,,Zj)dzj. 



Equation ( 14 ) giving the conditional distribution of the hitting scenario becomes 

Ylfll /{z,<y,c} ^(S.,M)(yr,, Z,)c?Zi 



The conditional distribution of the extremal functions Qt{y, t, ■) in Equation (15) is based 
on the kernel Pt(y, df). Using the existence of a Radon-Nykodym derivative for the finite 
dimensional margins of /i, we obtain 

y,/ s G dz) = dz. 

^t(y) 

This approach is exploited in ^ for Brown-Resnick max-stable processes. Indeed, the 
model turns out to be regular. 



5 Proofs 

5.1 Proof of Theorem [2] and Proposition [l] 

For the proof of Theorem [2| we need the following lemma giving a useful characterization 
of the i^'-extremal random point measure. If Mi, M2 G Mp(Co) are such that M2 — Mi G 
Mp(Co), we call Mi a sub-point measure of M2. 

Lemma 1. The K -extremal point measure $^ is the unique sub-point measure $ 0/ $ 
such that 

$ G and $ - <l G C^,(max(l>)). 

Proof of Lemma [l| First the condition $ — $ G C^(max(<l>)) implies 

max($ — $) <K max($) and max($) =k max($). 

Let / G [$ - I*]. The condition $ - $ G C;^(max(<l)) implies / <k max(<l). Since $ 
is a sub-point measure of max($) < max($) so that / <k max($) and / is i^T-sub- 
extremal in $. 

Conversely for / G [$] , the condition $ G implies the existence of to ^ such that 
/(to) = max($)(to)- Hence /(to) = max($)(to) and / is i^-extremal in $. 

□ 
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Proof of Theorem [2] : First note that $J(Co) = if and only if $ = 0. This occurs 
with probabihty exp[— /^(Co)] and in this case $J = = 0. The first claim follows. 
Next, let k > 1. According to Lemma [l| $J(Co) = if and only if there exists a /c-uplet 
{4>i, . . . , 0fc) e such that 

k k 

eC+ and $-5^5,^, eQ(vti0,)- 

1=1 i=l 

When this holds, the fc-uplet (0i, . . . , 0^) is unique up to a permutation of the coordinates 
and we have 

k k 

^5^,=$+ and $-^5^^ = $^. 

i=l 1=1 

Hence the sum 

„ fc-1 

•JU j=i 

is equal to ,$^,)gAx_B} ^tci^o) — ^ otherwise. Using this and Slyvniak's 

formula (see Appendix A.l), we get 

fc!P[($+,<l>^^)GAx5, $+(Co) = fc] 

fc-1 

This proves Theorem |2] □ 

Proof of Proposition [ij In the case yu(Co) < +oo, $ and a fortiori $J are a.s. finite. 
Suppose now yu(Co) = +oo, so that <l> is a.s. infinite. If inft<^K vi't) = 0, then there is 
to & K such that ri{to) = (recall rj is continuous and K compact). This implies that 
0(to) = for all G [$] and hence = $ is infinite. If inftg/^ r]{t) = e > 0, then the 
support of $J is included in the set {/ G Co; f{K) > e}. From the definition of Mp(Co), 
this set contains only a finite number of atoms of $ so that $ J must be finite. □ 

5.2 Proof of Propositions [2] and [3] 

Proof of Proposition [2| According to equation ([T]), for all x > 0, 

P[r^(t) = x]= exp[-/2t(x+)] - exp[-/it(x)]. 
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and f[i](t) = 0] = exp[— /it(0+)]. The equivalence between (ii) and (iii) follows. 

The equivalence between (i) and (ii) is a consequence of Theorem [2] with K = {t}, k = 1 

and A = B = Mp(Co): we get 

J Co 

exp[-fit{y)]fit{dy). 



'[0,+oo) 

It remains to prove that this probability is equal to 1 if and only if (ii) is satisfied. To 
this aim, we compute 

'm(Co) = 1] = / e-'^li^yp.^y^ydxntidy) (20) 

^(0,+oo)2 



/ e ''fit{A^)dx, 

J(0,+oo) 



'(0,+oo) 

where A^ = {y > : Jitiy) < x}. Since Jit is cag-lad, non-increasing and tends to 00 at 
0, Ax = (inf Aj., cxd) 7^ for all a; > 0. Furthermore using equation (20) and the fact 
that fit{Ax) < X, we get that P[<l>|'^|(Co) = 1] = 1 if and only if fitiA^) = x for all x > 0. 
We see easily that this is equivalent to condition (ii) and this completes the equivalence 
between (i) and (ii). 

We now prove Equation (Isl). Assuming that conditions (i)-(iii) are met, it holds 



i+ G A) = P[<|.+j e A $+}(Co) = l] 
with A = {6f, f G A}. Theorem |2| with K = {t}, k = 1 and B = Mp(Co) entails 

e^) = £l{5,eAnc+j}P[*nC^-^(/)]/x(rf/) 
l{/eA} exp[-/ii(/(t))] fM{df). 
This proves Equation (|8|). □ 



Proof of Proposition [3[ First note that the inequalities (10) characterize the hitting 
scenario. Let r = (ri, . . . , r^) G Vk and define the sets 

a = {(/i,...,/£)ec vjGii,^!, f,>r^ V/i'}- 

and 

e 

Cr = {Y,6f^eAucoy, (/i,...,/,)Ga}. 



Note that Cr C and that 6 = r if and only if $J G C^. 
Furthermore, B = r and {<ff, ■ ■ ■ , v?/) G A if and only if G Aj. with 



Ar = {Yl e ^p(Co); (/i, . . . , /^) G a n a}. 
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Hence the following events are equal 

{6 = r, ...,vi)e A, <l>^. eB} = {$+ e A,, <l>^ e B, $+ (Co) = i} 
and Theorem |2] implies 

P[e = r, G A 

= I l\^^^^s,^^^^/[^eBnC-^{yUf,)] i,^\df„...,dfe). (21) 
• Jcf, ^ 



Finally, Yl^j=i^f3 ^ if and only if there exists a permutation a of such that 

(/o-(i), . . . , /(t(^)) G a n Cr- Such a permutation is unique and this proves the equivalence 



of Equations (11) and (21). □ 



5.3 Proofs of Theorems |3|, |4| and |5] 

Proof of Theorems [s] and [4| Note that r](t) can be expressed in terms of the hitting 
scenario and the extremal function as follows. For r G Vk, define the mapping F^- : Cg — )■ 

[0, +00)'= by 

Tr (/i, ...Ji) = (yi, ■■■,yk) with yi = fjiU) if U e r,-. 

Definition (|2| entails that for all t G 9j, r](t) = ^'j{t). This can be rewritten as ?7(t) = 
Fe(0]^, . . . , 0/). Using this, the probability 

P(r, A, 5, C) = P[e = r, (y.+, ...,¥.+) G A, G B, r/(t) G C] 

can be computed thanks to Proposition [3j 

P{r,A,B,C) 

= P[0 = r, (y.+, . . . , ^;) G A n F;1(C), G fi, r/(t) G C] 

= l^^^,^^^ HvMm V,,. /.}^ [^^Bnc, {vUm /.«^(d/„ . . . , rf/,). 

Now for each j G |1, i], we condition the measure f^{dfj) with respect to /^(t^.): Equation 



(12) entails 

P(r,A5,C) (22) 

= X /. Hv^-i^'^^^ /^>. v.. ^ 5 n (vj.,/,)] 

®-=i{A./yr,,f^/i)/it./f^yr,)} 
= / / iweM, /,(M)<y,c}P[$e5nC-(y)] 

®-=i{A./yr,,f^/i)/t^t./f^yr,)}. 

In the last equality, we use the fact that fjitrj) = yr^ a.s. under Pt^^{yT-.,dfj), whence 

{vj G liJl f, >r, yy^.fr] = {vj G [1,^1, /, v,.^,/,.} 

= {VjGll,£l, /,(t..)<y,c 
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andC^(V^ti/,) = C7,-(y). 

We now prove Theorem [s} Setting A = Cq and B = Mp(Co) in Equation (22), we 
obtain 

r, r^(t) G C] 

Using the fact that P[$ G C^(y)] = exp[— /it(y)] and performing integration with respect 
to ®j=iA^. (Ytj, c(fj), we obtain Equation (13) and this proves Theorem 3 




We now consider Theorem 4 Combining Equations (13)-(16) together with Equation 



(22), we get 




r, (^+,...,^+)gA $;,g5, v{t)eC] 
<cJA P[$GC,-(y)] ^t(y,r,rf/,...d/.K(rfy) 

= 11 Rt{y,r,{fj),B)Qt{y,T,dfi - ■ ■dfi)'Kt{y,T)ut{dy). 

Jc J A 

In particular, with A = Cq and B = Mp(Co), we obtain the relation 



(23) 



r, r^(t) G C] 



7rt(y,^)t't(c?y) 



c 



characterizing the regular conditional distribution 



r I ?7(t) = y] (see Appendix A. 2) 



This proves that Equation (14) provides the regular conditional distribution P[0 = r | 
ri{t) = y]. Similarly, Equation (23) entails that the regular conditional distributions 
m^t) ^ ■ I ^(t) =y,'S> = r] and G ■ | r]{t) = y, 6 = r, (</.+) = (/,)] are given 

respectively by Qt(y, t, ■) in Equation (15) and -Rt(y, t, ifj), ■) in Equation (16). 



We briefly comment on these formulas. The fact that the distribution Qtiy,T, 



in Equation (15) factorizes into a tensorial product means that the extremal functions 
iff, . . . ,(p^ are independent conditionally on ^^(t) = y and B = r. The fact that the 



distribution i?t(y,T, (/j), ■) in Equation (16) does not depend on r and (fj) means that 



conditionally on ?7(t) = y, is independent of 9 and {(pi, . . . , ^^^(q-))- The distribution 
i?t(y, ■) can be seen as the distribution of the Poisson point measure $ conditioned to lie 
in C^iy), i.e., to have no atom in {/ G Co; fit) <^ y}. It is equal to the distribution of 
a Poisson point measure with intensity l|j(t)<y}Ai('^/)- D 

Proof of Theorem [5} Remark that 

Ms)<z} = {(<|.+ ,<I.^OeC-(z)xC-(z)} 

= U {e = r,<^+(s)<z,...,(^+(s)<z,$^, gC-(z)} 
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where (z) is defined in Theorem El Using this, Theorem H entails 



F[r]{s) < z I r]{t) = y] 

J2 P[0 = r, ^+(s) < z, . . . , ifiis) < z, G C-(z) I r/(t) = y] 

reVK 

E ^t(y,r)gt(y,r,{/(s) < z}Oi?t(y,C-(z)). 



reVK 

The result follows since 



Q.(y, r, {/(s) < z} ) - n P.,^(y,,.{/(t,.)<y,.}) 



and 



GC-(z)nC,-(y)] 



exp[-/x({/(s) ^ z or /(t) ^ y})] 

exp[-/i({/(t) ^ y})] 
exp[-M{/(s) ^ z, /(t)<y})]. 



□ 



5.4 Proof of Propositions [4] and [5] 

Proof of Proposition |4f This is a straightforward application of Theorem |4] and [5j 
Take into account that when K = {t}, Vk is reduced to a unique partition of size C = 1 
so that = {t} and y?^ = 4>t- D 

Proof of Proposition [sj According to Proposition |4| P[0^ G ■ | T]{t) = y] is equal 
toPt{y, ■). For any measurable A C Cq and -B C (0, +oo), we compute 




1 



B J Co 



J Co 

oo 



1 



f(t) 



feA}'^{rgmB}fit)^idf)r ^dra{dg) 



Co ^0 



Co 



l{rfeA}l{rf{t)eB}r ^dra{df) 
^{feA}l{f{t)eB} Kdf) 



The second equality follows from the change of variable r = rg{t)/f{t) together with the 
relation g{t)a{dg) = 1. This proves that 



Ptiy,A)= I l{^f^Ayfit)a{df). 

Jco 
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According to Equation (17) 



P[77(s) < z I vit) =y] = Ptiy, {/(s) < z}) exp[-/i({/(s) ^ z, fit) < y})]. 



We have 



P,(y,{/(s)<z}) 



and 



/i({/(s) ^ z, /(t) < y}) 



This proves Equation (19). 



Co 



l{^fis)<.}f{t)a{df) 



l/rf^s^y, r.fi'^w„}r '^dra{df) 



JCo 

oo 



Co 



V 



{rf{s)^z,rf{t)<y}l 

dra{df) 



^co ^mif^^'^m 



i=l 



y 



-) <df). 



□ 



A Auxiliary results 
A.l Slyvniak's formula 

Palm Theory deals with conditional distribution for point processes. We recall here one 
of the most famous formula of Palm theory, known as Slyvniak's Theorem. This will be 
the main tool in our computations. For a general reference on Poisson point processes, 
Palm theory and their applications, the reader is invited to refer to the monograph [15] 
by Stoyan, Kendall and Mecke. 

Let Mp(Co) be the set of locally-finite point measures on Cq endowed with the 
cr-algebra generated by the family of mappings 

^^ N{A), A C Co Borel set. 

Theorem 6 (Slyvniak's Formula). 

Let ^ be a Poisson point process on Co with intensity measure fi. For any measurable 
function F : eg X Mp(Co) [0, +oo), 



E 



k-l 



L ) • • • ) 'TKJ 



[ E[F(A,...,/,,$)]/x^'=(#i,...,rf/,). 
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A. 2 Regular conditional distribution 



We recall here briefly the notion of regular conditional probability (see e.g. Proposition 
A1.5.III in Daley and Vere- Jones [3j). Let {y,Q) be a complete separable metric space 
with its associated a-algebra of Borel sets, {X, F) an arbitrary measurable space, and vr a 
probability measure on the product space {X xy,J^®Q)- Let denote the Af-marginal 
of TT, i.e. 'Kxi^) = 7r(v4 x 3^) for any A G J-". Then there exists a family of kernels K{x, B) 
such that 

- K{x, ■) is a probability measure on (3^, Q) for any fixed x E X\ 

- K{-,B) is an J-'-measurable function on X for each fixed B E Q; 
{Ax B) = J^K{x,B)7ix{dx) for any A e and B eQ. 



- TX 



These three properties define the notion of regular conditional probability. When tt is the 
joint distribution of the random variable (X, F), we may write K{x, ■) = ¥{Y & ■\X = x). 

The existence of the regular conditional probability relies on the assumption that 
3^ is a complete and separable metric space. Furthermore, for any J-" ® ^-measurable 
non-negative function / on X x F, it follows that 



f{x,y)TT{dx,dy) = / / f{x,y)K{x,dy)TTxidx). 
xxy J X Jy 

The following Lemma states the existence of the kernel 

{Ps(x,rf/); xG [0,+oo)'\{0}} 



satisfying Equation (12). This is not straightforward since the measure fj, is not a prob- 



ability measure and may be infinite. 

Lemma 2. The regular version of the conditional measure fi{df ) with respect to /(s) G 
[0, +oo)' \ {0} exists. It is denoted by {Ps(x, (i/); x G [0, +oo)' \ {0}} and satisfies 



Equation (12). 



Proof : Let | ■ | denote a norm on [0, +oo)'. Define A = {/ G Co; /(s) ^ 0} and, for 
i>0,Ai = {fe Co; {i + 1)"^ < |/(s)| < i'^} with the convention 0"^ = +oo. Clearly, 
A is equal to the disjoint union of the AiS. We note //*(■) = /x(- fl Ai) the measure on the 
complete and separable space Co U {0}. Equation (Q ensures that /Xj is a finite measure 
(and hence a probability measure up to a normalization constant) and there exists a 
regular conditional probability kernel Pg(x, rf/) with respect to /(s) = x. We obtain, for 
all F : [0, +oo)' x Co, 

/ F(/(s),/)/i(rf/)= / / F(x,/)i^(x,ci/)/.,((ix), 

where = {x G [0,+oo)'^; (i + 1)^^ < |x| < i^^}. Let us define Ps(x, c?/) a probability 
measure on Co by 

p,{^,df) = Y,h^^A,}P:i^,df). 

i>l 
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If F vanishes on {0} x Co, we obtain 



F{f{s),f)fi{df) 



J2 [ F(/(s),/)/i(rf/) 

J2 i f i^(x,/)p:(x,rf/)/i,(rfx) 

i>0 "^^^0 

^ F(x,/)P,(x,d/)/i,(dx). 



[0,+oo)*\{0} ^Co 



This proves Equation (12). □ 



A. 3 Measurability properties 

Lemma 3. $^ and $^ are measurable from (fi, J-", P) to (Mp(Co), A^p). 

Proof : From Definition [T| it is enough to prove that that the events {0, <k i]} G J-' 
and {(pi ftx f]} J-'-measurable. Let Kq be a dense countable subset of K and note 
that (f) <K f] if and only if there is some rational £ > so that < ri[t) — e for all 
t e Kq. Hence, for all n G N U {+oo}, 

<K V \ N = n} = y Pi {TV = n ; < T^{t) - e} 

= u n = ^ ' '^^w < '^^■w - ^> (24) 

e>Ot&Ko j<n 

and {0i <A' v} = U^o{*^* ^ ' N = n} E J^. Note the union over e is countable since 
e is taken rational. □ 

Lemma 4. The set and C]^{g) are measurable in (Mp(Co), A^p). 

Proof : Let g he a, continuous function defined at least on K and consider the Borel 
set 

A = {/ G Co; / g} c Co. 
The set C]^{g) defined by Equation ^ is equal to 

C^ig) = {Me Mp(Co); M{A) = 0} 

and is A^p-measurable. 

In order to prove the measurability of defined by Equation ([s]), we introduce a mea- 
surable enumeration of the atoms of a point measure M (see Lemma 9.1. XIII in Daley 
and Vere- Jones ^). One can construct measurable applications 

K : Mp(Co) ^ N U {oo} and : Mp(Co) ^ Co, i> 1, 

such that 

k(M) 

M=Y1 ^^'(M), M G Mp(Co). 

i=l 
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A point measure M does not lie in if and only if it has a if-subsextremal atom. 
Hence, 

+00 k 

Mp{Co) \C+=[j [j{K{M) = k- ^M) <K max(M)}. 

k=Oi=l 

Similar computations as in Equation p4| entail 

{k{M) = k- ^i{M) <K max(M)} = |J fl U^'^ = ^ ' ^^(^) < " 

whence is A^p-measurable. □ 
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